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ABSTRACT 

The recent ROSAT X-ray detections of hot intergalactic gas in three groups of galaxies are reviewed 
and the resulting baryonic fraction in these groups is reevaluated. We show that the baryonic fraction 
obtained, assuming hydrostatic equilibrium, should depend, perhaps sensitively, on the radius out to which 
the X-rays are detected, and the temperature profile of the gas. We find that the NGC 2300 group has 
a baryonic fraction out to 25' of at least 20%, thus over five times higher than in the original analysis of 
Mulchaey et al. (1993), and also much higher than one would obtain from big-bang nucleosynthesis, but 
similar to the other two groups as well as rich clusters. With this baryonic fraction, groups would be fair 
tracers of the distribution of baryons in the Universe if Sl/i| = 0.3. A baryonic fraction that increases 
with radius is consistent with the X-ray data from all three groups. However, a detailed analysis of the 
NGC 2300 group shows that the dependence of baryonic fraction on radius is not well constrained by the 
data, in part because of uncertainties in the estimated background. 

Subject headings: galaxies: clusters of — galaxies: intergalactic medium — cosmology: observations — 



1. INTRODUCTION 



While clusters of galaxies commonly emit diffuse X-rays 
associated with a hot intergalactic medium, as evidenced 
numerous X-ray observations with the EINSTEIN Ob- 
rvatory (e.g., Forman & Jones 1982), the presence of de- 
ctable diffuse hot X-ray emitting gas in small groups of 
laxies does not seem as common. EINSTEIN observa- 
tions have revealed diffuse hot gas in four compact groups 
(Biermann, Kronberg & Madore 1982; Biermann & Kro- 
*s£berg 1983; Bahcall, Harris k Rood 1984), but failed to 
^Ttetect X-rays in 2 other compact groups, while a third 
Compact group emitted X-rays that were probably asso- 
ciated with the individual galaxies (Bahcall et al.). Re- 
' ce^nt higher-sensitivity ROSAT observations with the Po- 
C9t ion Sensitive Proportional Counter (PSPC) have con- 
i^Kmed that some groups of galaxies emit X-rays originat- 
ing from a hot diffuse intergalactic medium. Indeed, short 
flS^ lOmin) snapshots from the ROSAT survey mode has 
Produced 12 X-ray detections among Hickson's (1982) 100 
jeompact groups and these detected X-rays are probably of 
■watergalactic origin in roughly 9 of these 12 compact groups 
Xfl. Ebeling, private communication). Deeper ROSAT 
PSPC images were recently obtained in four groups: the 
loose groups NGC 2300 (Mulchaey et al. 1993, hereafter 
MDMB), NGC 5044 (David et al. 1993), and the Hickson 
compact groups HCG 12 (H. Ebeling, private communi- 
cation) and HCG 62 (Ponman & Bertram 1993). In only 
in one group (HCG 12) was there no intergalactic X-ray 
emission. The analysis of the X-ray properties of the NGC 
2300 group has revealed a relatively low fraction of mass 
in gas and stars (hereafter, baryonic fraction). Mulchaey 
et a/.'s analysis produced a baryonic fraction of 4%, and 
their upper limit was 15%. On the other hand, Ponman 
& Bertram's (1993) baryonic fraction was > 13%, while 
David et al. (1993) get 10%, both within 250 kpc (through- 
out this Letter we assume Ho = 50 km s _1 Mpc -1 ). 

The baryonic fraction in clusters has been used as a 
constraint on the cosmological density parameter £1: if 



clusters are representative of the whole Universe, one has 
£1 = S!j/ fi, (White 1992). Standard big-bang nucleosyn- 
thesis produces a baryonic fraction on the scale of the 
whole Universe of 0.04 < tt b hl < 0.08 (Steigman 1989), 
where h^o = Ho/(b0 km s _1 Mpc -1 ). The baryonic frac- 
tion in clusters of galaxies is typically around 25% (e.g., 
White 1992, and references therein). One thus obtains 
SI ~ 4Slj ~ 0.24. As noted by MDMB, the implications of 
a baryonic fraction as small as S!j is that if the Universe 
has a density near closure (Slo — 1), groups of galaxies 
would be fair tracers of the baryonic content of the Uni- 
verse. 

Another way to set constraints on SI is to assume that 
groups are fair tracers of the ratio of mass to blue light 
in the Universe. The luminosity density of the Universe 
is such (Loveday et al. 1992) that the closure mass-to- 
light ratio is M / L = 780 /150, and in general one can write 
SI = M/L/(780h 50 ). The typical M/L of groups (e.g., 
Gourgoulhon et al. 1992), computed with the assumption 
that groups are in virial equilibrium, yields SI = 0.07, while 
a value of SI = 0.3 is obtained with the correction for 
the non-virialized cosmo-dynamical state of these systems 
(Mamon 1994). 

In this Letter, we consider how the baryonic fraction in 
groups (and clusters) of galaxies should vary with radius, 
assuming simple models for the properties of the hot X- 
ray emitting diffuse gas. We analyze the published data 
for NGC 2300 using a similar model to that applied to rich 
clusters to reevaluate its baryonic fraction at the present 
limiting detection radius, and we indicate which future X- 
ray observations could place sufficient constraints on this 
parameter. 

2. baryonic fraction versus radius 

The X-ray surface brightness profile of groups and clus- 
ters are usually well-fitted by a law of the form 



S(R) = S [l + (R 2 /Rl)} 



2M 1/2-3/3 



(1) 
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produces a central peak in the surface brightness profile. 
Here [3 is an empirical rather than a physical parameter 

Inverting the Abel equation relating X-ray surface 
brightness to the 3D emissivity in the spectral passband of 
the instrument, noting that the emissivity of a hot plasma 
varies as n 2 A(T), where A is the cooling function, call- 
ing n = dlnA/dlnT, and assuming an equation of state 
T ~ n 7_1 (i.e., isothermal for 7 = 1 and polytropic oth- 
erwise), one obtains (see Cowie, Henriksen & Mushotzky 
1987) 
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T(r) 



n(r) = n [l + (r 2 /Rl)} 

T [l + (r 2 /^)]"' (7 " 1) 
3/3/2 



(2a) 
(2b) 

(2c) 



l + n( 7 -l)/2 
From equation (2a), the total mass in gas can be written 

M gas (r) = Airn R 3 c iJm p M(x) , (3a) 
x = r/R c , (36) 
y 2 dy 



M(x) 



(3c) 



/o (l + y 2 ) 

where /j, is the mean particle weight in units of the proton 
mass m p . 

1 



(x 2 + 1) 



M(x) : 
M(x) : 

M(x) = sinh _1 a; ■ 



1/2 



-sinh 



tan 



for 6 
for 6 
for 6 



1/2 ,(4a) 
1, (4b) 
3/2 .(4c) 



(x 2 + 1) 1/2 

Writing the equation of hydrostatic equilibrium, the to- 
tal mass (that binds the gas) is (Cowie, Henriksen & 
Mushotzky 1987) 

kTr /dlnn dlnT N 

Mtot : 



2j6 



Gfj,m p 
kT R. 



!lnr 



!lnr 



(5) 



Giim p (1 + x 2 y+ s (-v-i) ' 

Using equations (3) and (5), the gas fraction can be writ- 
ten 

_ 2*G» 2 mjn Rl | ^ l+6H _ r) M(x) ^ 

3 ybkTo ' ' x 3 ' 

Figure 1 shows the gas fraction versus radius for isother- 
mal models with different 6s, scaled to T = IkeV, R c = 
100 kpc, and no = 10 -3 cm -3 . In general, the gas fraction 
should be multiplied by 

(no/10 -3 cm- 3 )(i? c /100kpc) 2 (T /lkeV)- 1 . (7) 

From Figure 1, one sees that for 6 = 1/2, the gas frac- 
tion increases sharply with radius, while for 6 = 3/2, the 
reverse is true, and finally for 6=1 the gas fraction is 
asymptotically constant, and equal to 2.3 times its core 
radius value. From equation (6), the gas fraction varies 
asymptotically as 

f g ~ (r/R c ) 2 -^ 2 -^ (8) 

so that a constant asymptotic gas fraction implies 6 = 1 for 
isothermal models and 6 = 3/2 for j = 5/3 polytropes. 
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FlG. 1. Scaled (see eq. [7]) gas fraction (eq. [6]) vs. radius 

for isothermal models. 



Figure 2 shows the same as Figure 1, but for poly- 
tropic models (all with 6=1). Again the gas fraction 
should in general be multiplied as in equation (7). For 
non-isothermal gas distributions, the gas fraction increases 
with radius outside of the central region. 
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FlG. 2. Scaled (see eq. [7]) gas fraction (eq. [6]) vs. radius 

for 6 = 1 polytropic models. 



Note that certain parameters (6 = 1/2 with j = 1 or 
6 = 1 with 7 = 5/3) yield unphysical gas fractions over 
unity at large radii (R > 100i? c ). 
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3. COMPARISON WITH OBSERVATIONS 

The ROSAT X-ray observations are not strongly con- 
straining for NGC 2300. For example, in Figure 3, we 
plot the intergalactic X-ray surface brightness profile of 
the NGC 2300 group. The points are from MDMB (where 
we omit the three points at 18', 21' and 24' for which 
there is obstruction from the window support structure) 
and uncorrected for the background while the curves are 
fits to these points using the model in equation (1), with 
a constant background Sb g added everywhere. Of course, 
MDMB have determined the background independently of 
their NGC 2300 observations, but their adopted value of 
Sbg = 7.0 x 10 -4 arcmin _2 s _1 (D. Burstein, private com- 
munication), could be uncertain by 10% or more, and this 
is why we allow it to be a free parameter. Figure 3 clearly 
shows that very different good fits (without the point at 
39', which is well above the surrounding points and may 
be contaminated by a point source) provide very different 
6 and core radii. 
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FlG . 3 . Surface brightness profile of the NGC 2300 group in- 
tergalactic X-ray emission. Data (uncorrected for background) from 
Mulchaey et al. (1993) are shown as points, while the best fits 
to isothermal-/? models (eq. [1] plus a constant background) shown 
in the top three lines of Table 1 are respectively shown as dotted, 
solid, and dashed curves. Points at 18', 21', and 24' are omitted 
as they suffer from obstruction from the window support structure 
ofthePSPC. 

Using the results from §11, the top three rows in Table 
1 illustrate how these three fits return different total mass 
(M t ot, eq. [5]), gas mass (M gas , eqs. [3] and [4]), gas frac- 
tion (f g , eqs. [4] and [6]), and baryonic fraction (/j) within 
25' of the center of the diffuse gas in the NGC 2300 group. 
The fourth row in Table 1 lists what would be obtained if 
we adopted the MDMB parameters of Sbg = 7.0 x 10 -4 , 
f3 = 0.85 (J. Mulchaey, private communication), and in 
subsequent rows, we investigate alternative fits imposing 
either R c or [3 in addition to Sb g . As MDMB, we as- 
sume isothermality and we convert gas fractions to bary- 
onic fractions using M ga i(25') = 6 x lO 11 /ig O 1 M . The gas 
number density necessary to normalize the gas fractions 
(eq. [3a]) is obtained by writing the bolometric luminosity 
as 



L M = ln 2 A(T)dV=A7rn 2 R 3 c A(T ) 



-Rmax/-Rc 



x 2 /(l+x 2 ) 2S dx 



where A is the Raymond and Smith (1977) cooling func- 
tion. The density calculation uses the cooling func- 
tion evaluated for the emission weighted temperature 
(0.86 keV) and is 8.6 x 10 -24 ergs cm 3 s _1 (using a heavy 
element abundance of 6%, the best fit value of MDMB) 
and i? max = 25' is the extent of the X-ray emitting hot 
gas. 



TABLE 1 

Reanalysis of the NGC 2300 Group for R < 25' 
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Notes: (1) Background in 



10" 

4 . 



arcmin 

-2 , 



(2) Cen- 



tral surface brightness in 10 arcmin s . (3) Core ra- 
dius in arcmin. (4) Surface brightness shape (eq. [1]). (5) 
Goodness of fit. (6) Bolometric luminosity in 10 42 ergs s _1 . 
(7) Central gas density in 10 -3 cm -3 . (8) Gas mass within 
25' in 1O 12 M . (9) Total mass within 25' in 1O 12 M . (10) 
Gas fraction within 25'. (11) Baryonic fraction within 25'. 



Table 1 shows that the best fits (rows 1 to 3) conspire to 
a surprisingly robust gas mass within 25'. Moreover, the 
fit with the MDMB parameters is significantly worse than 
the three fits shown in Figure 3, because their adopted 
background is too high. The best fit baryonic fractions 
are > 20%, over 5 times larger than the value quoted by 
MDMB (4%). In fact, our analysis shows that even the 
MDMB parameters should return a baryonic fraction as 
high as 11%. The discrepancy with MDMB's low result is 
due in part to MDMB probably computing the baryonic 
fraction within 15' (as stated in the main body of their 
paper) while in their Note added in proof, they mention a 
total mass one-third lower within a region extending two- 
thirds further (25'). For comparison, NGC 5044 has a 
baryonic fraction of ~ 10% within 250 kpc (David et al. 
1993) and HCG 62 > 13% within a similar region. Also, 
David et al. detect diffuse X-ray emission out to 400 kpc. 
Since NGC 2300 fills the PSPC field of view, the extent of 
its gas is > 0.5 Mpc. 

The first two groups have low values of 6: Ponman's & 
Bertram's analysis of HCG62 yields 6 = 0.54, while David 
et al. find 6 = 0.79 for NGC 5044. Also, the tempera- 
ture profiles for NGC 5044 and HCG 62 are slightly non- 
isothermal outside of the cooling radius with effective j of 
1.1 and 1.15 respectively. These values of 5 and j imply 
that the gas fraction in HCG62 and NGC 5044 is increas- 
ing with radius near the limit of the ROSAT detections. 
On the other hand, the NGC 2300 data is not constrain- 
ing enough to evaluate 6 (from eq. [2c], 6 = 3/3/2 for 
the assumed isothermal model). Indeed, decent fits yield 
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8 = 0.65 to 8 > 2.5. The high background used by MDMB 
forces a high asymptotic slope for the surface brightness 
profile, yielding values of 8 much higher than that found 
for the other groups and for clusters. 

4. DISCUSSION 

Groups are probably more representative of the Universe 
than are clusters, as they include altogether perhaps five 
times as many galaxies as do clusters. With a baryonic 
fraction over 20%, the inner 25' (330 h^kpc) of the NGC 
2300 group would be a fair tracer of the baryonic fraction 
in the Universe if £1 ~ 5 £lj ~ 0.3/ig 2 . These numbers 
are similar to the constraints from group M / Ls, with (if 
Ho = 50) and without (if Ho = 100) corrections for non- 
virialized states, respectively (see §l). 

From equation (8), a constant baryonic fraction can be 
reached with 8=1 isothermal gas or with 8 = 3/2 poly- 
tropic gas with index j = 5/3. Figure 1 shows that 
for isothermal 8=1 gas, the asymptotic gas fraction is 
roughly 2 to 2.5 times larger than at 25'. Therefore, an 
extrapolation of the gas fraction to large scales will yield 
a baryonic fraction /j > 40%, and will thus be consistent 
with baryonic nucleosynthesis if SI < 0.15 /ig 2 . Taken at 
face value, the observations of the three groups discussed 
here suggest 8 < 1 and j > 1, both of which imply bary- 
onic fractions that increase with radius, which extrapo- 
lates at large scales to even lower values of SI. However, 
as noted in §111, high values of 8 can produce equally good 
fits to the NGC 2300 data, and thus produce constant or 
decreasing baryonic fractions. Note that a similar trend of 
increasing gas fraction is found in clusters (e.g., Durret et 
al. 1993). 

To reconcile the data for groups with SI = 1, as favored 
by inflation, one requires a baryonic fraction that decreases 
at large radii to an asymptotic value near S!j ~ 0.06 /ig 2 , 
hence implying 8 > 1 (isothermal gas) or 8 > 3/2 (7 = 5/3 
polytropic gas). Although these values of 8 seem inconsis- 
tent with the X-ray observational data of the three groups 
studied here, there are two possibilities to remain consis- 
tent with inflation: 1) While the gas fraction is nearly con- 
stant within the group, it could fall to zero outside of the 
group, so that groups of galaxies would be sites of higher 
baryonic fraction in the Universe, i.e. biased tracers of 



the distribution of the baryonic fraction in the Universe; 
2) The cosmological constant, A (not to confuse with the 
cooling function used above with the same symbol) is non- 
zero, i.e., SI = 0.2 and A = 0.8, which is consistent with 
the distribution of cluster temperatures deduced from X- 
ray observations (Bartlett & Silk 1993). 

The model surface brightness profiles of the NGC 2300 
group have a similar shape in the inner region, and are 
thus indistinguishable, although uncertainties in the PSPC 
background subtraction allow these profiles to diverge sig- 
nificantly beyond R > 45', as seen in Figure 3. This un- 
certain behavior at large radii for NGC 2300 and the in- 
creasing baryonic fraction in HCG62 and NGC 5044 make 
clear that it is of fundamental importance to observe X- 
ray emission at relatively large distances from the centers 
of groups to determine the extent of the gas. 

As a final note, if the three groups studied here are still 
in the stages of cosmological collapse (as is argued for all 
non-compact groups by Diaferio et al. 1993 and Mamon 
1994), then the equation of hydrostatic equilibrium used 
here (eq. [5]) may not apply, for two reasons: 1) The gas 
may not have time to react to the rapid changes of the 
global potential of the group; 2) Even if the gas follows 
the potential of the whole group, one needs to add a term 
d(pv r )/dt in equation (5). Ponman & Betram (1993) con- 
clude that the compact group HCG 62 is past full collapse 
and currently in the process of slowly shrinking by orbital 
energy dissipation via dynamical friction. In general, it is 
reasonable to assume that during the collapse of a group, 
the denser inner parts will collapse earlier and be close to 
equilibrium today, the question then being how far out is 
the gas in equilibrium. It may well be that if hot inter- 
galactic gas found in a group presents unusual properties, 
this could reflect a departure from hydrostatic equilibrium 
possibly caused by the group's cosmological collapse. 
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